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Abstract 

The vibration of the dimer excited by STM current on Si(OOl) surface is investigated. We 
describe this system by the Hamiltonian which has the electron-vibration couphng term as 
the key ingredient. In order to characterize the transition rates induced by STM current be- 
tween vibrational states we have introduced the effective temperature of the vibration which 
differs from the temperature of the substrate. The behavior of the effective temperature 
depends on the substrate temperature and STM current in highly nonlinear manner and 
qualitatively changes around 50K of the substrate temperature. At lower temperatures, the 
effective temperature strongly deviates from the substrate temperature and reaches a few 
hundred Kelvin for the typical values of STM current. At higher substrate temperatures, the 
effective temperature reduces to the substrate temperature. On the basis of these behaviors 
of the effective temperature, we solve the puzzle of the symmetric-asymmetric crossover in 
dimer images of STM observation in the ordered state of c(4 x 2). 
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1. Introduction 

Intensive studies have been made on the atomic structure of the reconstructed Si(OOl) 
surface both experimentally and theoretically. On Si(OOl) surface, neighboring atoms form 
buckled dimers which have two stable tilting angles. Along the (110) direction, the dimer 
rows is formed by the dimers. The principal features of the surface reconstruction on Ge(OOl) 
surface are similar to Si(OOl) surface. The ground state of the structure of the dimer arrange- 
ment of both surfaces is c(4 x 2) and the ordered c(4 x 2) phase turns into a disordered state 
at room temperature through the order-disorder phase transition. ^"^^^ In the ordered state, 
the dimers are observed as the asymmetric images by scanning tunneling microscopy (STM). 
In the disordered state, the symmetric appearing images are observed by STM, and these 
symmetric appearing images have been attributed to the rapid repeat of the orientational 
change (flip-flop motion). 

In the STM observation, semiconductor substrates are doped into n-type ones or p-type 
ones for the electric conductivity. On Si(OOl) surface, serious differences in the property of 
the dimer system between the n-type substrates and p-type ones have not been reported at 
temperatures higher than 60K. 

Recently, p(2 x 2) phase on Si(OOl) surface-*^^"^^^ is observed on highly doped n-type 
substrates below 40K by STM. A phase transition is also observed by low-energy electron 
diffraction^^) (LEED) at 40K on highly doped n-type substrates. On the p-type substrates, 
c(4 X 2) phase on Si(OOl) surface is observed down to 9K.^^) 

At much low temperatures, the symmetric STM images of p(2 x 1) structure^^'^^'^^^ are 
observed at some conditions of the tip bias voltages Vt and the tip currents It in STM 
observations. Yokoyama^^^ reported that STM images of the dimers at 5K on boron doped 
(B-doped) p-type Si(OOl) substrate at the condition of Vt = ±1V and It = 50pA appear to 
be symmetric, although most of the dimer images on the substrate are well asymmetric to 
make the ordered c(4 x 2) structure at 63K, and that the symmetric STM images at 5K are 
observed also on Antimony doped n-type Si(OOl) substrate. Kondo^^^ reported subsequently 
that the STM images at 20K on phosphorus doped n-type Si(OOl) substrate at the condition 
oiVt — 2.0V and -1.6V, and It — 300pA appear to be symmetric. These observations of the 
symmetric images in the ordered phase are puzzling and have not been explained yet. 

The first-principles calculations (FPC) with high accuracy^'''''^^) show that the ground 
state on Si(OOl) surface is c(4 x 2). The results of FPC show that the energy of the elec- 
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tronic system decreases, and the energy of the lattice system increases, as the tilting angle 
of the dimer increases. The equilibrium angle is determined by the competition between 
the energies of two systems. Thermal equilibrium properties on the surface at temperatures 
higher than lOOK have been successfully explained by Monte Carlo simulations (MCS) with 
Metropolis algorithm based on the results of FPC^'^^ by mapping the two stable orientations 
of the dimer onto the Ising system. The time-resolving dynamical Monte Carlo sim- 
ulations (TDMCS) performed for the dimer systems of Si(OOl)^^'^^'' surface and 00(001)^^-* 
surface show that the most of the long range fluctuations of the dimer arrangement on Si (001) 
and Ge(OOl) surfaces are carried by the phase boundary of the c(4 x 2) domain in the dimer 
row at low temperatures. In order to perform TDMCS on Si(OOl) surface and on Ge(OOl) 
surface, the model potentials for the continuous values of the tilting angle are obtained based 
on the FPC rcsults^^'^*^) . The TDMCS results with the model potentials reproduce well the 
results of the time-resolving STM obscrvations^^'^^^. Since at temperatures higher than 60K, 
as mentioned above, serious dopant effect are not reported on Si(OOl) surface, the reported 
MCS results based on the results of FPC remain valid. 

The surface localized states of 7r-band and 7r*-band^'^'^^'^^'^^'^°~^^^ are formed mainly 
through the hybridization of the dangling bonds on the dimers. FPC^'^'^^^ show that the 
energies of the top of the vr-band and the bottom of 7r*-band stay a few ten meV above 
from the top of the valence band of Si and about SOOmeV below from the bottom of the 
conduction band of Si, respectively. The band widths are obtained to be about leV in both 
of TT-band and 7r*-band. These features are observed by angle-resolved UPS^°^ (ARUPS) and 
scanning tunneling spectroscopy i^'^^'^^) (^gxs). On the n-type substrates, the Fermi level of 
the system stays near the bottom of the conduction band, which is much higher than the 
bottom of 7r*-band. The large band bending appears near the surface region of the n-type 
substrate at much low temperatures. On the p-typc substrates, the Fermi level stays near 
the top of the valence band. Therefore, the top of the 7r-band is near the Fermi level. The 
large band bending does not appear near the surface region on the p-type substrates. Be- 
cause the dopant effect on the substrate is not taken into account in the FPC, the results of 
FPC reproduce well the surface properties of p-type substrate, but some modifications are 
necessary for n-typc substrate at much low temperatures. 

In the present study, we take the scope to the vibration excited by STM current on the 
p-type substrate. We will introduce the Hamiltonian which has the electron-vibration cou- 
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pling term. The excitation rates by the STM currents for the vibration will be obtained. The 
rates will be characterized by introducing the effective temperature of the vibration. The 
mechanism for the symmetric- asymmetric crossover in dimer images of STM observation on 
Si (001) surface will be presented. 

2. Model 

We introduce the Hamiltonian H representing the coupling between the electronic system 
and the localized vibrational system, 



H — He + He-Y, 

He = E£fe4cfe + 5^£pcJcp + ri(X]7fc4)($Z7pCp) + H.c. I 

k P { k P J 

= X]^°4ca + Saa^'a + (Er^a'^Ca + H.C.) + Y,^p4>'^p+ (X^Tpa^Cp + H.c), 

a a P P 

i/e-v = huj{b^b +7:)+ Se{b^ + b)a^a, 

^ ] \^k\ ^k-i 
k 

k 

Tp = r7p, (1) 



where Cfc, Cp and b are the annihilation operators for the electronic states in the vr-band, 
for the states in the conduction band of the STM tip, and for the vibrational state in the 
rocking mode of the dimer on Si(OOl) surface beneath the STM tip. We assume that the 
vibrational states are the harmonic oscillator, and the energy of the state with the vibrational 
number n is liuj{n-\- 1/2). The energy of the vibrational state Tiuj is measure by the electron 
energy loss spectroscopy^^^ as 20meV. o) = ■jkcl and J2p 7pc], are the creation operator 
for the spatially localized states \a) in the surface dimer and in the prominence of the 
tip, respectively. The tunneling current is through the the spatial localized states. The 
coefficients of the spatial localized states are normalized as J2k l7feP — 1 X^p ITpP — 1- 
The diagonalized term of J2k^kCk^k is rewritten by a and Cq as the first three terms in 
the third line of eq. (1), where Cq is the annihilation operators for the electronic states 
orthogonalized with \a). The amplitude of the vibration modifies the one electron energy 
of I a), and H has the electron- vibration coupling term, because the energy of the electronic 
system, as mentioned above, decreases as the tilting angle of the dimer increases. The term 
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of 5e represents the electron-vibration coupling. 5e is the couphng constant between the 
amphtude of the vibration and the energy of the one electron state of \a). The vibration 
is excited or deexcited through the electron- vibration couphng by the STM current. The 
similar Hamiltonian of eq. (1) has been used in the precedent studies of vibrations of 
molecules adsorbed on metal surfaces^^"^^^ . In the present study, however, the Hamiltonian 
is applied to totally different system, the surface localized semiconductor state on Si(OOl) 
surface. 

In the present study, we treat the case that the tip bias voltages 14 is positive and so 
large that the Fermi level in the conduction band of the STM tip is below the bottom of the 
TT-band of the p-type Si(OOl) substrate; T4 is assumed to be shghtly larger than IV, because 
the band width of the 7r-band is about leV. The elastic inter-band transition of electrons 

from the vr-band to the conduction band of the tip is induced by the presence of the tip. 
The elastic inter-band transition rate dcias of electrons from the surface to the tip without 
coupling to the vibration^"^""^^^ is given by the T-matrix element and the Fermi distribution 
functions: 



27]- 

C^elas = 2— ^ |r*G'a(£)rc,|^ 

^exp(/3(£«-£F)) + l^^^ " exp(/3(£, + eyt-£F)) + l^^^^" ~ 

^ ^exp(/3(£„ -£f)) + 1 l + exp(-/3(£p + eyt - £f)) 

Ga{e) = (a|(£ + iO+-//e)"'|a> 

1 



£-£a-A(£)+iA(£)' 
A(£) = A,(£)+At(£), 

As(£) = 7r5^|r,,|25(£-£«), 

a 

He) - r de'^, (2, 



TT 

where £f is the Fermi level of the electronic system, Ga{£) is the Green function for the 
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electron system, A{s) is the width of the projected density of states for |a), As{e) and At{e) 
are the components of A{e) in the surface and in the tip, respectively, P denotes the Cauchy 
principal value, and E^^ and D^^ are the energy level of the top of the 7r-band and the width 
of the TT-band, respectively. In the present study, — (-Btt — ^f) is assumed to be of the order 
of a few ten meV and D^^ is assumed to be large enough PD^ 3> 1. The temperature of the 
system is assumed to be so low that the broadness of the Fermi distribution around ep in 
the surface state and around ep — eVt in the tip state arc not contributed to the evaluation 
of the integration in eq.(2). We assume that At(£a)/As(£a) is small in the STM observation 
in the typical conditions enough to be approximated as A(£a) As(£a)- Because in the STS 
experiment at 5.5K, the tunneling spectra of the 7r-band are observed to be almost symmetric 
band, we assume that stay at the center of the 7r-band in energy, and A{ea) = 0. Therefore, 
(Teias is approximately obtained as 



(Xeias ~ —At{ea)A;\ea)D^. (3) 

The inter-band transition of electrons from the 7r-band to the conduction band of the tip 
induced by the presence of the tip can excite the vibrational state in the dimer through the 
electron-vibration coupling. The inelastic inter-band transition rate (J^qX^ of electrons^^"^^^ 
from the surface to the tip coupling with the excitation from the ground state of the vibra- 
tional number n = to the excited state of n = 1 is obtained in the lowest order in \Ss\ 
as 



= 2^j:\¥,MHe-.W,0)ME{p')-hu;-E{a')) 

a' ,p' 



-Dn+E„ 



^1 1 

^ ^exp{f3{£a - £f)) + 1 exp(/3(£p - hu + eVt - £f)) + 1 

r Ps{ea)Pt{ea)L)n, 



a' 

1 , 



7r'{e-ea-He)y + A^{e)) 
Pt(£) = ^|(p'|a)|^(5(£-£p/) 
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= 1( ^ ) (4) 

where E{a') and E{p') are the one electron energies of \a') and \p'), respectively, and ps 
and pt are the projected density of states for \a). \a') and \p') are the stationary states of 
He connected to jo;) and \p), respectively. Namely, the stationary states are given as jo;') = 
|a) + {£a + iO+ — He)''^ra\a) and \p') — \p) + (e^ + iO+ — He)~^Tp\a). The projected density 
of states are approximated in the same way as eq. (3): Psisa) ~ (7i"As(£a))^^ and ptisa) ~ 
At(£a)(7rAg(£a))~^ The inelastic transition rate ctq'^^ is given by use of the approximation 
for the projected density of states as 

cr'oii ^ ^\{ea)A:%Sa)D,. (5) 

The elastic transition rate (Teias and the inelastic transition rate ctq^Xi depend critically 
on It and do not on T, because in eqs. (3) and (5), At depends critically on the spatial 
tunneling gap of STM. The fraction of the inelastic transition rate to the elastic transition 
rate^*^"^^^ W — cr^^i/aeias, however, does not depend on At, and is obtained from eqs. (3) 
and (5) as 

^ - '^^ 

Eq. (6) shows that the fraction of cr^t^Xi/crcUs is given dominantly from the intrinsic property 
of Si(OOl) surface, and independent of the STM current. In the present study, as mentioned 
already, Vt is assumed to be so large that the Fermi level in the conduction band of the tip 
is below the bottom of the 7r-band, and hcu <S eVt. Prom these assumption, the inelastic 
inter-band transition rate ct^'^q of electrons from the surface to the tip coupling with the 
deexcitation from the vibrational number n = 1 to n = is well approximated to take the 
same value as o"q°^]^. Thus the inelastic inter-band transition rate c^^ln-ti of electron from 
the surface to the tip coupling with the vibrational transition from n to n ± 1 is given as 
<1n±i = niax(n, n ± l)a^Xi- 

At finite temperatures T, the deexcitation of the vibrational number n = 1 to n = 
through the inner-band excitation of the electron-hole pair creation in the 7r-band can occur. 
In the inner-band excitation, both the initial state and the final state of the electron are 
within the states of |q;'), which have dominant component in the 7r-band on the surface. The 
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deexcitation of the vibration with the inner-band excitation docs not occur at OK. The rate 
of the deexcitaion ctJ^^q of vibration through the inner excitation is independent of It and 
depends on T: 



2tt 

= 2—Y.\KAH.-.\a',,l)\\5{E{a'^)-hw-E{a[)) 



— / de ps{e)ps{e + hu}) 

1 1 

^ ^exp(/3(£-£F)) + l ~ exp(/3(£ + flu; - £f)) + 1 



TT 1 — exp(— 



/3;i^n + exp(/3(E,-;ia;-£F))^' 
The transition processes of cXeias, cri'^o; '^o^^i; '^I'^o; "^o^i schematically shown 
in Fig. 1. The essential framework of the formalism for cXeias, Cq^d and ct"^q is sim- 




Figure 1: Fig. 1 The transition processes of the inter- 
band transitions and the inner-band transitions. The 
rates of the inter-band transitions of CTeias, <^T^qi and 
o'o'^i between the 7r-band and the conduction band of 
the STM tip depend on It, and not on T. The inner- 
band transitions of <t™o and CTg"^]^ within the 7r-band 
denpend on T, and not on 1^ . 



ilar to that of theoretical studies^^ *^^^ on the STM induced vibration of molecules ad- 
sorbed on metal surfaces. In the problems of molecules adsorbed on metal surfaces, tran- 
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sitions corresponding to the inner-band transitions in the present study have also been 
treated. The transitions corresponding to the inner-band transitions, however, do not 
exhibit any temperature dependence in the case of metal surfaces. Eq. (7) shows that 
(J^o depend strongly on T; nfJha'^^o ^ AWex.p{f3{E^ - Ef)) for -P{E^ - £f) > 1 and 
Trfdha'^^Q pa 2W for —(^{E^, — hu — Sp) ^ 1. The rate of the excitation a}^^^ of vibration 
through the inner-band deexcitation is easily obtained by the essentially same way as (jf^Q. 
a^i ^ iA/n)W[FiT)/iexp{phu;) - 1)]. 

The total transition rate of the excitation Ro-*i of the vibrational number n = to 1 and 
the deexcitation i?i_>o of the vibrational number n = 1 to are given by the sum of rates as 



D f ^ I ^inn I _.int 

exp[phLj) — 1 



^ -){To' + -WFiT))+Wa,i^, (8) 



'l-exp(-/3M ^ 
where tq is the life time of the vibration through the background phonon system at OK. 

The probability P„ to find the vibrational state at n satisfies the rate equation represented 
by Ro^i and -Ri_o: 



d-Pn 

= {n + l)i?i^o-Pn+i + ni?o-*i-Pn-i - {nRi-,0 + {n + l)Ro-,i)Pn, I < n, 

at 
dPo 

= Ri^oPi — Rq^iPq- (9) 

In the stationary distribution of P„, P„ is derived as P„ = r"(l — r) from eq. (9), where 
r = Ro^i/Ri^Q. 

The tip current It are given by the sum of the elastic inter-band transitions and the 
inelastic inter-band transitions: 



e 



n=0 n=0 n=l 

1 + r 

[1 + W- ]c7elas. 



(10) 
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When T, W, and It are given, i?o^i and Ri^o are obtained selfconsistently with eqs. (8) 
and (10). The distribution of the dimer vibration in the STM observation is characterized 
by the effective temperature Tef of the dimer vibration defined as 



r = exp(-^u;/(A;Bref)). (11) 

Throughout this paper, T denotes the temperature of the substrate. Tgf differs from T 
in general and depends on T and It nonlincarly. The key quantity for T^f is the fraction 
W — cr^Q^i/creias, which is derived from |fe| and Ag as shown in eq. (6). The couphng 
constant 5e between the amphtude of the vibration and the energy of \a), is approximately 
estimated as 



= , (12) 

oO e=eo 

where Oq (= 19.1° = 0.333rad ) is the equilibrium tilting angle of dimer in the c(4 x 2) 
structure, and is the mean value of the amplitude of the vibrational number n = 1. A9 
is obtained by lu and the effective value of the moment of inertia / as A^ = {3h/ (2/a;))^/^. As 
mentioned already, the model potentials for the continuous values of the tilting angles of the 
dimers on Si(OOl) surface are obtained based on the FPC results. We can get approximate 
values of / and dsa/ d9 at 9 = 9q by the model potentials: 

o9 $=$0 

' = ^- (13) 

where k is the value of the second partial derivative of the model potentials for the tilt- 
ing angle of the dimer beneath the tip at 6' = ^o- We derive the value of k from the 
model potentials, as A; = 1.1 x lO^meV/rad^. Using eq. (13), we get approximate values as 
7 = 1.2 X 10-23meVsVrad^ dea/d9 = -3.7 x lO^meV/rad, and A^ = 5.2 x IQ-^rad. Prom 
eq. (12), \6e\ is approximately estimated to be 1.9 x lO^meV. Prom the half width of half 
maximum of STS for 7r-band of Si(OOl) at 5.5K, Ag is roughly estimated to be 2.5 x lO^meV. 
Prom these values, the key quantity W is obtained to be about 0.6. 



3. Results and Discussions 
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In the present study, we take the scope on the p-type substrate. The Fermi level of the 
p-type substrate is well approximated to be at the center of the top of the valence band 
of the substrate and the acceptor level of dopant at low temperatures. At temperatures 
low but higher than a few Kelvin, is practically fixed at the top of the valence band. 
In the present study, E^^ — ep is approximated to be — 22.5meV for B-doped substrate and 
— 32.5meV for Ga-doped substrate, respectively. The life time of the vibration tq through the 
background phonon system at OK is expected to be of the order of 10~^s for semiconductors. 
We assumed for tq to be 8ns. 

The excitation rate i?o-»i and the deexcitation rate i?i_»o for B-doped substrate are shown 
in Fig. 2. Both rates are monotonically increasing function of T and It in highly nonlinear 
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Figure 2: Fig. 2 The total rates for the excitation of the dimer vibration and the deexcitation for 
B-doped substrate. The rates are obtained selfconsistently with eqs. (8) and (9). (a): the excitation rate 
Ro^i- (b): the dexcitation rate iii^o- 



manner. The rates scarcely depend on T for T < 20K, because the inter-band transitions 
are dominant in them and the inner-band transitions scarcely contribute to them. When 
T increases around 20K, the excitation and deexcitaion of the vibration through the inner- 
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band transitions of electrons become effective and the rates increase steeply for T. The 
deexcitation rate -Ri_>o increases much steeper for T than -Ro-^i around 20K ~ 40K. The 
rates scarcely depend on for T > 70K, because the inner-band transitions are dominant 
in them and the inter-band transitions scarcely contribute to them. 

Since the magnitudes of the rates are smaller enough than the frequency of the vibration 
4.8THz derived from hu = 20meV, the effective temperature of the dimer vibration Tcf = 
—huj/{kBlog{Ro^i/Ri^o)) well characterizes the distribution of the vibration. Tcf for B- 
doped substrate are shown in Fig. 3. At low temperatures of T < 20K, Tef scarcely depending 




Figure 3: Fig. 3 The effective temperature of the dimer vibration (Tef) 
for B-doped substrate. Tef reaches about 500K at It = 200pA at low tem- 
perature T < 20K. Tef decreases steeply around T f« 20 K ~ AOK, as T 
increases. 



on T, increases steeply with It, becomes about 250K at It — 50pA, and reaches about 500K 
at It — 200pA. The rapid flip-flop is expected in such high Tef. These high Tef appearing 
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at low temperatures of T < 20K, essentially explains why the symmetric STM images are 
observed at low temperatures. When T increases around 20K, T^i starts to decrease steeply, 
because of steep increase of -Ri-^o- ^ef becomes 70K around T Ri 40K and It ^ 150pA. At 
higher temperatures T > 60K, T^i increases gradually with T, scarcely depends on 1^ within 
the typical values of It, and takes the nearly same value as T. These low T^f explains why 
the asymmetric STM images recover in temperature range higher than 50K. 

Tcf for Ga-doped substrate arc shown in Fig. 4. The essential feature of Tcf shown in 
Fig. 4 is same as in Fig. 3. The steep decrease of Tef in Fig. 4, however, is starting around 
T 30K which is higher than for B-doped substrate shown in Fig. 3. The temperature 
where the inner-band transitions exceed the inter-band transitions in magnitude on Ga-doped 
substrate, is higher than that on B-doped substrate and the ratio of these two temperatures 
is almost the same value as the ratio of ep — E^^, 32.5/22.5=1.44, since cJ^^q, as mentioned 
already, is approximated ~ 4W cxp{j3{ET^ — Ep)) / [tt f3h) for — — sp) ^ 1- 

At low temperature T < 20K, the atom at the up-position beneath the tip is rapidly 
transferred to the down-position through the dimer flip-flop motion by quasi-thermal ex- 
citation of the vibration of high Tef. After the flip-flop motion to the down-position, Tef 
becomes much lower than at the up-position, because tip current is expected to be much 
smaller at down-position than at up-position. Because of the low effective temperature Tcf, 
the transition rate of the atom from the down-position to the up-position is much smaller 
than that of the opposite transition; the fraction of the staying time at the down-position 
is much larger than that at the up-position. The spatial STM gap between the prominence 
of the tip and the surface dimer is tuned for It to be the STM condition. The tip current 
It is time averaged one, and is tuned almost at the down-positon at low temperatures. In 
the scanning of STM tip at much low temperatures, the atoms beneath the tip are almost 
always tied to the down-positions, and the atoms at the up-position scarcely contribute to 
the STM images. As the result of the large fraction at the down-position, the symmetric 
dimer images of the p(2 x 1) structure is observed even in the ordered state of c(4 x 2). 
These feature that the tip ties the atom to the down-positions in STM observation on Si(OOl) 
surface corresponds to the experimental results reported by Mitsui and Takayanagi^^^ . 

In the present study, we demonstrate that the inner-band transitions play the impor- 
tant role to decrease Tcf at the temperatures of a few ten Kelvin where deexcitation of the 
vibration through the inner-band excitations of the electron-hole creation occur. This de- 
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Figure 4: Fig. 4 The effective temperature of the dimer vibration (Tef) 
for Ga-doped substrate. Tef reaches about 500K at It = 200pA at low 
temperature T < 30K. Tef decreases steeply around T « 30 K ~ 60K, 
as T increases. The decreases on Ga-doped substrate appears in higher 
temperatures than on B-doped substrate. 

excitation through the inner-band excitation can occur also when enough amount of holes 
in the 7r-band or enough amount of electrons in the 7r*-band is induced. For example, when 
the electrons and holes are induced by infrared photons or by other methods, Tef will be 
expected to decrease in much low temperature in STM observation. Here we propose the 
notion of "the excitation-induced coohng" , namely, the decrease of Tef essentially induced 
by the electronic excitation in vr-band or 7r*-band. 

4. Conclusion 

In order to explain the symmetric- asymmetric crossover in dimer image of STM observa- 
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tion on Si(OOl) surface, wc calculate the effective temperature Tef of dimer vibration. The 
crossover behavior is semiquantitatively understood by the nonhnear dependence of Tef on 
the substrate temperature and STM current. The symmetric image in the ordered state of 
c(4 X 2) results from the large fraction of the atom at the down-position induced by STM 
current. We have proposed the notion of "the excitation-induced cooling" under more gen- 
eral situations than that induced by the substrate temperature. 
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